Description of self-synchronization effects in distributed Josephson junction arrays 

using harmonic analysis and power balance 
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Power generation and synchronisation in Josephson junction arrays have attracted attention for 
a long time. This stems from fundamental interest in nonlinear coupled systems as well as from 
potential in practical applications. In this paper we study the case of an array of junctions coupled 
^~>' to a distributed transmission line either driven by an external microwave or in a self-oscillating 

mode. We simplify the theoretical treatment in terms of harmonic analysis and power balance. We 
apply the model to explain the large operation margins of SNS- and SINIS-junction arrays. We show 
the validity of the approach by comparing with experiments and simulations with self-oscillating 
es-SIS junction arrays. 



Josephson junctions (JJs) are well known to be able to generate microwave power. Their dynamics can also phase- 
lock either to an external microwave bias or to the self-generated signal from an array of junctions Q. A common 



O ' configuration is to couple a number of JJs through a distributed transmission line resonator ^5]. Both overdamped 

I , [^] and underdamped Q junctions in this configuration arc known to be able to mutually synchronize. Long arrays 

Oh' of JJs coupled to a transmission line driven by external microwave bias are known to exhibit properties, which can 

^ , only be explained by the collective dynamics of the J J-transmission line system 0, 0| . 

^ ; • Coupled Josephson junction arrays have been modelled by semiclassical dynamical simulations o^' 

, terms of lengthy perturbation analysis . The dynamics of coupled arrays is complex making engineering of such 

■ devices challenging, often based on trial and error. However, in most practical situations it suffices to study the 
I , situation, where all the junctions are phase- locked to a driving signal (either external or self-generated). Recognizing 

^ ' this, we develop a simplified modelling technique. We have previously considered the case of external microwave bias, 

, for which we calculated the gain in the limit of very small power self-generation Here we extend the analysis to the 

^ ' case of arbitrarily large gain, and apply it to self-stabilized overdamped JJ arrays. We analyze also the self-oscillating 

' , array by means of a power balance equation. Finally, we compare the results of the analytic model with simulations 

^ ■ and experiments. 

, We assume that the tunnel element is effectively voltage biased. This is ensured by a small impednace |^(/)| ^ 

^\ ^of /Ic, connected in parallel with the tunnel element (see Fig. 1(a)). Here Ic is the critical current, $o is the fiux 

C , quantum and / — lu/2tt is the drive frequency. We further assume that the dynamics is phase- locked to the driving 

\l • signal. Below we drop the explicit /-dependence from the function notation, since all the quantities are evaluated at 

I that frequency. The condition for phase-locking (to the first Shapiro step) is 

o 

where ii = |Z|/i/$o/, ^^-^o ~ {h — ^of/R) and Ji is the Bessel function of the first kind. Here 7^ is the DC bias 

I ' current, R is the DC resistance in parallel to the junction and Ii is the amplitude of the driving signal. 

Using the assumptions above, we get the fundamental frequency component of the Josephson current Iji (see Fig. 
1(a)) Within harmonic approximation the voltage across the tunnel element is Vji (Ji) = Z (/i — Iji (/i)), where 

Q , the phase is referred to the driving signal h. We define dynamical impedance Zd (/i) = Rd + iXd ~ V.ji (/i) //i, in 
series with ujLt ^ X^, where Lt is the transmission line inductance per junction length. We neglect Xd, assuming 

• I— I . Xd < ujLt. The effects discussed here follow from the real part i?d (-^i) = Re (Zd (/i)) = Re (Z) — Re (.Z/ji//i). Re (Z) 

■ is the real part of the linear impedance parallel to the junction, and Re(Z/ji/Ji) is due to the Josephson effect. 

^ , Knowing Rd, various properties of arrays are readily obtained. We consider first the case of external microwave 

■ - - ' drive (Fig. 1(b)). With \Rd\ -C ujLt the (nonlinear) attenuation per junction is 



\SIo\ < Ic 



Ji (n) I , (1) 



n I ^ Rd {h {xk)) 
ilKhKXk))^-- , (2) 

where Zq is the transmission line impedance. If Rd > 0, the junction attenuates and if Rd < 0, it amplifies the 
driving current. This is a generalization from the theory of low-loss transmission lines (e.g. Il2j V The current 
distribution Ii (xk) along an array {xk being the position of the fc:th junction) can be recursively obtained from 
Ii (xfc_|_i) = Ii (xk) exp (?7 (/i (xk))). This also generalizes our previous results _j5] to arbitrarily large gains. 
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FIG. 1: (a) A microwave current biased junction, (b) A Josephson array coupled to a matched transmission line with external 
RF source, (c) A Josephson array coupled to a A/2-resonator. An nX/2 resonator is a straightforward generalization, (d) JJ 
with resistive environment, (e) JJ with capacitive environment. 



Our second case is a resonator driven by self-generated power of the JJ:s (Fig. 1(c)). We note that the generated 
(or dissipated) power of a JJ at position is —Rd {Ii (a;fc)) Ii {xk)'^ /2. The current distribution in the system can 
therefore be solved by means of a power balance equation 

1 ^ 

-^Y.^d{Il{xk))xh{xkf =Pe.,, (3) 
k=l 

where the left side is the generation due to all N JJ:s and Pex is the excess loss, i.e. all power loss mechanisms not 
related to the junction or the shunt. In a high-Q resonator with \Rd\ , \Xd\ -C ujLt the resonant frequency determines 
the generated frequency. The position dependence of Ii (x) depends also only on the resonator properties, whereas 
the overall amplitude depends on the gain and dissipation of JJ:s and on Pex- The criterion needs to be satisfied 
for all Ji (xk)- 

We identify two limits for Rd in terms of junction type. The first is a resistive environment (Fig. 1(d)): Z — R 
at the drive frequency. This is relevant to SNS or SINIS junctions with appropriate parameters. In this case we get, 
using above definition and results of Ref. 's^ for Rd,R ^ Rd{Z = R) 



in which J{ is the derivative of the Bessel function. 

The second limit is a capacitive environment, Z « 1/jujC, which is valid for unshunted SIS junctions. The subgap 
leakage resistance can typically be neglected. This is also a good approximation for externally shunted SIS (es-SIS) 
junctions, if \R + ujLs\ ^ 1/ujC (Fig. 1(e)). Here Lg is the series inductance of the shunt resistor. We get 



i?d,c(/i) = Re(Z)-^-^, (5) 

where Re(Z) « for unshunted junctions. For es-SIS junctions the losses in the shunts can be the dominant mechanism 
of linear dissipation, although their effect on JJ oscillation (second term in ||SJ)) is small. For the arrangement of Fig. 
1(e), Re (Z) « (Leff/Lsf i? if > L^ff, where L^ff ^ Lt - l/u'^C. 

It is experimentally known that long overdamped arrays driven by external microwave compensate for the trans- 
mission line attenuation Q . This is the reason why SNS and SINIS Josephson voltage arrays work in spite of the huge 
attenuation due to intrinsic damping 0,^3- With the above notation, self-compensation occurs if the amplitude of 
the propagating signal is initially between Ji^min and Ii^max such that (see (0, 0) 

r,(/i) = 0,/i 

riih) > 0,/i,,„i„ < /i < /{, (6) 

.max ■ 
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FIG. 2: Operating margins of an overdamped J J array in self-stabilizing mode according to Q and Open circles (left 
axis): Maximum step width as a function of '^of /IcR- Gray-shaded area (right axis): pump current range maximising the step 
width. 



In addition, /i.min E^nd /i^max need to satisfy (Q. Then Ii {x) approaches the stable point ![ towards the end of 
the transmission line. Fig. 2 shows the maximum step width 25Jo,max and pump current margin /i.max — ^i.min, 
within which the array self-stabilizes the microwave current. Typical values of the design parameter ^of/IcR of 
overdamped arrays vary between 0.5 and 2. In this range, the maximum self-stabilized step width varies between 
0.2/c and 0.6/c- Using the parameters of Fig. 2 in Q (<I>o///ci? — 1.3), our model predicts 2i5/o,max ~ 400 fj,A and 
201og]^Q (/i, max/A, min) ~ 18 dB, which are in excellent agreement with the exprimental data (2(5/o,max ~ 400 fiA and 

201ogio (/l,max//l!min) > 10 dB). 

In SIS and es-SIS junctions, rj does not change sign as function of Ii (see the second term in ((SJ). Therefore the 
tendency for power self-stabilization is not present. The much smaller attenuation of the capacitive environment 
makes phase-locking of all junctions possible in this case for a reasonable number of junctions ^fkil3|. 

Finally, to study self-oscillating arrays, we prepared arrays embedded in open-ended A/2 or A microstrip resonators. 
The circuits were fabricated using the Nb/Al/AlOx/Nb fabrication line of VTT li'l and measured at 4.2 K. The 
resonators were 800 fim or 1600 /im long Nb strips having a resonant frequency / nominally about 80 GHz. The 
resonators had 2-28 es-SIS JJ:s centered around the node(s) of the standing wave. A photograph of one of the 
structures is shown in Fig. 3(a). Measured IV curves are shown in Fig. 3(b). 

The AC current distribution in the resonators is given as Ii (x) = Jg sin (27ra;/A). The excess loss is given as 
Pex = loLtqIq/'^Q. where Lto is the total inductance of the resonator and Q is the intrinsic Q-value. Using Eqs. 
(EJ , lO the power balance stands ^ ^ 

f=i [{Leff/Lsf Rx {Iosm{2TTXk/X)f - 2$o/ x SIqJ = TrfLtoIo/Q- This is 
solved for /q. For 6Iq < 0, no physical solution exist, since Rd,c can be negative only for positive SIq. We get 



h [x] = . / 5 Tj sm — — . (7) 

V 27r/W2Q+(l/2)(£e///is)'i?Ef=iSin'(2^^feA) \ ^ J 



SIo < Ic 



Ji 



( h{xk) 



(8) 



Here we have rewritten ^ for convenience. For the resistive environment also the negative side of the step 51^ < 
appears. 

The solutions for which all junctions are phase-locked to the driving signal appear in the IV curves of Fig. 3 as 
steps with voltage (V) — N^of (for simplicity we omit the small dependence of / on the bias point due to Xd). This 
happens for the solutions of iQ that satisfy for all Xk- Solutions where no phase locking occurs are described by 
a resistive IV curve determined by the shunt resistors. We assume here that these are the only types of solutions. 

The theoretical IV curves obtained in this way are shown as solid lines in Fig. 3(b). The fitting parameters are 
the Q-value and the resonant frequency /. The Q-value was fitted to the case with N = 2, for which the effect of the 
parameter spread is the smallest. The result Q = 420 was used for all the arrays. This corresponds to a transmission 
line loss of a = tt/XQ x 8.7 dB « 41 dB/m, very close to that of an earlier measurement (50 dB/m) The loss 
is dominated by the dielectric (PECVD Si02) loss. The resonance frequencies for all devices were very close to the 
nominal value 80 GHz. Other parameters were obtained independently either from the IV curves or from the known 
circuit and material properties. 

The fact that steps appear only for SIq > as expected is clarified for = 2 in Fig. 3(b). The step widths are 
close to the values predicted by the theory, though smaller. The difference is explained by parameter spread (the 
shunt resistances varied up to ±5% within a chip, changing the effective bias point i5/o by tens of /zA). 
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FIG. 3: (a) An experimental array with a A/2-resonator having N = S JJ:s. Typical circuit parameters are Ic ~ 280 nA, 
C = 59 pF, i? = 0.26 Q., Ls = 3.4 pH, Lt = 0.2 pH, Lto = 5.7 pH and / = 80 GHz. (b) Current-voltage curves measured from 
the resonators with A'' JJs (small red circles). The device with A = 2 x 14 JJ:s is based on a A-resonator. Others are based on 
A/2 resonators. The voltage is scaled to that of a single junction and curves are lifted by 40 /iV with respect to the one below 
for clarity. The analytic curves are shown by solid black lines and the simulated points are denoted by large blue diamonds. 



To further study the vahdity of the model we performed time-domain simulations using Stewart-McCumber model 
for the junctions and a distributed model of a dissipative transmission line [l^ to describe the resonator sections. 
We found good agreement with the analytic model. The model step widths are well reproduced, though the small 
bending of the plateau due to finite Xd is not accounted for by the single-frequency model. As assumed, in all cases 
the numerical solutions were such that either all the junctions were synchronized or nonsynchronized. However, for 
simulated devices having even larger number of JJs per resonant node than our experimental devices, we found also 
partially synchronized states. Also in these cases the bias range of the completely synchronized state was correctly 
predicted by the analytic model. 

In summary, we have developed a simple analytic technique for quantitative modelling of distributed JJ systems with 
a strong preferred frequency set by either an external RF signal or a high-Q resonance. The results clarify differences 
between underdamped and overdamped systems, a source of some controversy recently. The model provides an 
intuitive way of understanding the dynamics of phase locking and power generation in JJ arrays. The model predicts 
results in agreement with simulations and experiments and it can be used to engineer practical devices, such as 
mm-wave local oscillators and JJ voltage standards. An interesting alternative voltage standards is a JJ array with 
no external microwave source, locking the self-generated signal directly to a frequency reference. The power balance 
equation is easily generalized to externally loaded systems to calculate the power output. Other types of physical 
systems can be analyzed similarly. For example, the exact dual for systems described here is a distributed parallel 
array of quantum phase slip junctions 0|, where the "conductor gain" of J J arrays is replaced by a "dielectric gain" . 
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of Excellence in Low Temperature Quantum Phenomena and Devices) and EU through RSFQubit (FP6-3749). The 
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